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We consider the recent RHIC data on the transverse single spin asymmetry (SSA) An, measured 
in p^p — > 7T° X processes at mid-rapidity by the PHENIX collaboration. The measurement is con- 
sistent with a vanishing SSA. We analyze this experimental information within a hard scattering 
approach based on a generalized QCD factorization scheme, with unintegrated, transverse momen- 
tum dependent (TMD), parton distribution and fragmentation functions. It turns out that, in the 
kinematical region of the data, only the gluon Sivers effect could give a large contribution to An; 
1 its vanishing value is thus an indication about the possible size of the gluon Sivers function (GSF). 

Approximate upper limits on its magnitude are derived. Additional constraints obtained combining 
available parameterizations of the quark Sivers function and the Burkardt sum rule (BSR) for the 
Sivers distributions are also discussed. 

in 



O 
O 
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> ! I. INTRODUCTION AND APPROACH 

i— I , 

■ We have recently discussed a hard scattering approach to hadronic interactions, based on the assumption of a gener- 
alized QCD factorization scheme which involves unintegrated TMD parton distribution and fragmentation functions; 
thepartonic intrinsic motions are also fully taken into account in the elementary perturbative QCD (pQCD) processes 
d S 9 • This scheme has been applied to the computation of unpolarized cross sections for single- inclusive particle 
production, pp — > hX, at high energy and (moderately) large pt [H, and to the computation of the transverse SSA 

cu: a _ d(? - do-± 

^ An - dcTTTdcTi (1) 

in p 1 p -> h X processes 0, H H H . 

Transverse single spin asymmetries can originate, even with a short distance helicity conserving pQCD dynamics, 
. \ from spin-fcj^ correlations in the soft components of the hadronic process, the distribution and fragmentation functions. 

There could be many such correlations. However, the study of SSA based on the assumption of a generalized 
factorization scheme 0, Q shows that the correct treatment of the elementary pQCD dynamics, with all phases 
related to the non-collinear and non-planar partonic configurations, leads to a strong suppression of all contributions, 
with the exception of the Sivers 0], and, to a lesser extent, the Collins [6] mechanisms. 

The issue of the validity of the factorization scheme with unintegrated partonic distribution and fragmentation 
functions is still an open one. Such a scheme has been shown to hold for Semi-Inclusive Deep Inelastic Scattering 
(SIDIS) and Drell-Yan processes @, BQ, while it is not yet clear whether or not it holds for inclusive one-particle 
production in pp processes; in such a case it is difficult to account for the gauge links necessary to ensure the gauge 
mvariance of the TMD distribution functions 0, El, E3, 0, El ■ At this stage we consider our factorized description 
oi pp — > 7r X processes as a phenomenological model based on a natural extension of the usual collinear approach 
for the same process, and of the factorized scheme with unintegrated partonic distributions proven for SIDIS and 
Drell-Yan processes. 

The Sivers function has recently received a lot of attention: data on azimuthal SSA in semi-inclusive deep inelastic 
scattering (SIDIS) processes from the HERMES collaboration at DESY [l5[ and from the COMPASS collaboration 
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at CERN lial have allowed, for the first time, a direct extraction of the Sivers functions for u and d quarks inside a 
proton 

mm, Hindus. simiiari y. 

u and d Sivers functions have been extracted from purely hadronic processes [l[ , 
in qualitative a gree ment with the SIDIS results. Also the possibility of accessing the gluon Sivers function has been 
investigated P.l22l. l23l l24l]. Although one might expect that spin-fc_i_ correlations are stronger for valence quarks - 
as the large xf data from the FNAL-E704 collaboration [25[ and the STAR collaboration at RHIC-BNL [26| seem 
to indicate - one knows that gluons play a dominant role in many high energy hadronic processes; it would be very 
interesting to see whether or not the gluon density inside a transversely polarized proton depends on the intrinsic 
motion. 

We address here the issue of the largely unknown gluon Sivers function, A N f g / p i (x, k±), and its possible contribution 
to the SSA Am for the p^p — > tt° X process, in the framework of the generalized factorization scheme, with pQCD 
elementary dynamics, of Refs. [H, H, Q- In this approach, the general structure of the cross section for the polarized 
hadronic process (A, Sa) + (B, Sb) — > C + X, is given by [|[ 

Ecda (A,S A ) + (B,S B )^C+X , dXadxhdz 

" C) 



^ / : ~ ~ 2 d2fc ^« d 2 k± b d 3 k± c 8(k±c ■ Pc) J{ki 

I _L07T X dXjy Z S 



a,b,c,d.{\} 



X P\ /A jf A fa/A,S A ( x a.,k_La) P\^)f B fb/B,S B ( x b,k± b ) (2) 

a 5 a b ' b 

x Mx c ,x d , K ,x b Mt, ctXd , K>K S(s + t + u) D^{z, k ±c ) , 

where all parton intrinsic motions are fully taken into account, both in the soft, non perturbative components and in 
the hard, pQCD interactions. 

In Ref. [3J] the structure of Eq. © was extensively discussed; its main features are the appearance of several spin 
and k± dependent distribution and fragmentation functions (with a partonic interpretation) and the non-collincar 
partonic configurations which lead to many k± dependent phases. As a consequence, it was explicitly shown how the 
integration over the parton intrinsic momenta leads to strong suppressions of most contributions to the unpolarized 
cross sections and to the transverse single spin asymmetry An. The only sizeable contributions to the latter, in the 
kinematical region (large positive xf) of the E704 [25[ and STAR data [20] come from the Sivers and, less importantly, 
from the Collins mechanisms. The dominance of the Sivers effect is even more pronounced in other Xp ranges, namely 
at xf < 0. 

In Ref. Q also the flavour decomposition of the Sivers effect was performed; while the quark contribution is 
totally dominant at large and positive xp values (for polarized protons moving along the positive Z-axis), the gluon 
contribution may be sizeable at xf ~ 0, the mid-rapidity region. This can be easily understood, since x™ ln (the lowest 
kinematically accessible value of the light-cone momentum fraction of parton a inside the transversely polarized initial 
proton) must be larger than xp\ An at large Xf is then mainly driven by valence quark properties. Indeed, an analysis 
of the E704 data allowed a first extraction of the Sivers functions for u and d quarks 1] . In principle, also inclusive 
production in the negative xp region might be sensitive to small x a gluons inside the polarized proton (being hit by 
large x b partons inside the unpolarized one). We shall further comment on this point after Eq. ([5]). 

Data in the mid-rapidity region are available from the E704 [27| and PHENIX [28| experiments. The kinematical 
regime corresponding to negative values of xp has been covered by the STAR collaboration (29| . Preliminary results 
for charged pions are also available from the RHIC-BRAHMS experiment [30| . 

In these kinematical regions A^(p^p — * n X) is largely dominated by the Sivers effect alone, and Eq. @ gives [l|: 



E ~W~~ %r~ - X! / dXadXbd 2 Z d 2 kj_ a d 2 k ±b d 3 k JL7r 5(kj_ n ■ p c ) J(fc ±7r ) 

" P-n ^ Pit abed ^ ^ a ^ 5 



da ab ^ cd 

x Af a/p r(x a , k ±a ) f b / p (x b , k ±b ) s 2 — (x a , x b} s,t,ii) S(s +t + u) D„/ C (z, k ±7T ) , (3) 

where (M denotes the proton mass) 

k± a 

A /a/pT (x a ,kx a ) = / a/pT (x a ,kj_ a ) - f a/pl (x a ,kj_ a ) = A N f a/p1 (x a ,k ±a ) cos0 a = -2 — j- f 1T (x a , k± a ) cos0 a . (4) 

A N fa/pi ( x a> k± a ) [or f^p(x a , k± a )] is referred to as the Sivers distribution function of parton a inside a transversely 
polarized proton. <j) a is the azimuthal angle of the intrinsic transverse momentum k± a of parton a. We follow here the 
notations and kinematical conventions of Ref. [H, with the polarized proton moving along the positive Z-&xis, in the 
pp cm. frame; the observed pion is produced in the XZ plane, with positive X values; spin f and J, are respectively 
along the positive and negative Y-eods. 
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Eq. ([3]) gives the numerator of An, Eq. (1); the denominator is just twice the unpolarized cross section, which is 
given by 

E^da unp ^ [ dx a dx b dz 2 2 3 „ 

— 75 ~ > / 5— a kj_ a d k±_ b d k± w 6(k± n ■ p c ) J(k± v ) 

d S p v ^ J TTX a X b Z 2 S 

a,b,c,d 

d(j a ^ — ' c< ^ 

x fa/ P (xa,k ±a ) f b / p (x b ,k ±b ) s 2 — — (x a ,x b ,s,t,u) 6(s + t + u) D^/^Zjk^) . (5) 

The azimuthal phase factor cos0 a appearing in the numerator of An, Eqs. ([3|) and (j4|), plays a crucial role and 
deserves a comment. The only other term depending on (j) a in Eq. ([3]) is the partonic cross section, in particular 
via the corresponding Mandelstam variable t. At large positive xf and moderately large px the (average) values 
of t are relatively small. Therefore, the (dominant) i-channel contributions, proportional to 1/t 2 , depend sizeably 
on <p a , so that the d 2 k± a cos <j> a /i 2 integration in Eq. ([3]) does not necessarily suppress An- Instead, for negative 
values of xp, all partonic Mandelstam variables are much less dependent on <f> a , so that one is roughly left with 
the d 2 k± a cos cf> a integration alone, which cancels the potentially large Sivers contribution. As a consequence, the 
possibility of gaining information on the gluon Sivers distribution from the recent STAR and BRAHMS data at 
negative values of xf is frustrated. Notice that, due to the much lower values of ^/s involved, this suppression caused 
by the cos0 a dependence would be much less effective for the kinematical regimes of E704 and of the proposed 
PAX [3l| experiments Q. However, E704 never measured An for negative xf values, while the PAX experiment is 
still in the proposal and planning stage. 

Similar considerations lead to a strong suppression of the gluon Sivers contribution to the SSA for p^p — -> 7 X 
processes at STAR, in the negative xf range and at a (pseudo) rapidity 77 ~ —4 [3^ |. 

The same arguments do not app ly to inclusive hadronic processes at mid-rapidity and moderately large pt values, 
for which data from PHENIX [28| are already available, for neutral pions and charged hadron production. For these 
processes, the gluon contribution is dominant and the Sivers effect can survive the phase integration. In the next 
section, we shall therefore consider in detail this case, aiming at a possible derivation, from data, of useful direct 
constraints on the gluon Sivers function. 



A. Constraints from the Burkardt sum rule 



Indirect constraints on the GSF could also be obtained from a sum rule for the Sivers distribution recently derived 
by Burkardt (33|. The BSR states that the total (integrated over x and k±) transverse momentum of all partons 
(quarks, antiquarks and gluons) in a transversely polarized proton must be zero, 

(k±) = J2 < fc i>" = f dx [ d 2 k^k 1 _Y, A A/ P T {x, k ± )=0. (6) 

a J J a 

Naively, at the partonic level considered in our approach, this sum rule simply corresponds to total (transverse) 
momentum conservation inside a transversely polarized proton, since the unintegrated distribution function 

fa/pi y^a: k± a ) fa/p O^as ^-_La ) + ^ A/a/pT ( Xa ,kj_ a ) , (7) 

introduced in Eq. has a clear probabilistic interpretation. On the other hand, it is well known [I3,[ll| that a non- 
vanishing Sivers effect requires initial/final state interactions, which might spoil the simple partonic interpretation. 
The Burkardt sum rule ensures the non-trivial result that momentum conservation holds also in this situation. The 
validity of the BSR has been explicitly verified in a diquark spectator model calculation of the Sivers distribution in 
Ref. [34[. In this paper, we consider the BSR as an additional theoretical tool in order to obtain indirect information 
on the gluon Sivers distribution, once the quark distributions are known. 

We notice, however, that a strict use of the BSR requires integration over the full x range, including the poorly 
known small x region, of the single Sivers functions, which might even result in divergences. Therefore, in the 
following we shall simply check whether or not the parameterizations of the Sivers functions jwhich will turn out from 
our phenomenological analysis of the data within the theoretical approach of Refs. [l], 0, [Hi fulfill the BSR in the 
limited x range considered, without extrapolating them to very low x values. 

In the next section we present our results on the GSF, based on the PHENIX [28[ experimental data on An(p^p — > 
ir° X). Further comments and conclusions are given in the last section. 
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II. PHENOMENOLOGY OF THE GLUON SIVERS FUNCTION 

We consider the PHENIX data [II] on A N for the p ] p — * tt° X process at RHIC, at ^fs = 200 GeV, with p T ranging 
from 1.0 to 5.0 GeV/c and mid-rapidity values, |r;| < 0.35. In this kinematical regime, at the lowest pt values, £™ m 
can be as small as 0.005. Therefore, partonic channels involving a gluon in the transversely polarized initial proton 
dominate over those involving a quark. This gluon dominance, together with the (almost) vanishing of all possible 
contributions to An other than the Sivers effect, allow to interpret the data - showing tiny values of An - in terms 
of useful constraints (upper bounds) on the gluon Sivers function, as it will be shown. As px grows, x™ ln increases 
and the dominance of the gluon channels becomes less prominent. 

Another set of data, for comparable rapidity and pr ranges, has been collected several years ago by the E704 
Collaboration |27j, at a lower energy, yfs ~ 20 GeV. In this case, however, even at the smallest pr values, x™ ln 
remains large enough so that gluon channels are not dominating; a possible mixing with quark initiated contributions 
then prevents to get clean constraints on the GSF from E704 lower energy data. 

Before analyzing in detail the reliable (for our purpose) PHENIX data, let us summarize what we know so far 
about the Sivers functions. As already discussed, the E704 and STAR collaborations have measured large SSA for the 
p^p — ► irX process at large positive xf and moderate pr values. In this kinematical regime, since x™ in > xf, gluon and 
sea-quark contributions should be negligible. Indeed, in Ref. [IJ we have shown that reasonable fits to the SSA can be 
obtained by using valence-like Sivers functions for u and d quarks, which turn out to have opposite signs. In addition, 
fits to the weighted azimuthal asymmetries measured for pion production in SIDIS with a transversely polarized target 
, Uj| , lead to independent comparable parameterizations of the u and d Sivers distributions [l7], 0, HU [Hj] . 

Notice that in all these phenomenological analyses the Sivers functions for gluons and sea-quarks have been assumed 
to be vanishing. That is a reasonable assumption, considering the kinematical region of the available data, and the 
lowest order decoupling of gluons in SIDIS. A natural question arising at this point is the following: how do the 
valence u and d Sivers functions alone, so far extracted, perform with the mid-rapidity PHENIX and E704 data? The 
answer is that they predict an almost vanishing SSA, compatible with both sets of data (see Fig. [TJ for the PHENIX 
results). Not only, but their parameterizations [l|,[l3] are also compatible with the Burkardt sum rule, Eq. ([5]), yielding 
(fc_i_) ~ within a 10% accuracy. It looks like there is no need to introduce other contributions to the Sivers effect in 
addition to that coming from valence quarks. 

Of course, this cannot be a definite, although simple and appealing, conclusion, as the data used for the extraction 
of the u and d quark Sivers functions are insensitive to the small x region. A large gluon Sivers function would not 
affect the analysis of the SIDIS, E704 and STAR data at large positive xf', however, it would strongly affect the 
description of the mid-rapidity PHENIX data. The real question is now: how much does the small value of A n 
measured by PHENIX suppress the gluon Sivers function? 

In order to answer this question we compute An according to Eqs. (UJ and imposing different conditions on 

the gluon Sivers functions, trying to understand what is the maximum value of \A N f g / p r (x, k±)\/2f g / p (x, k±) allowed 
by the PHENIX data. We follow Ref. [l(, adopting the same factorized gaussian k± dependence for distribution 
and fragmentation functions and the same parameterization for the Sivers distributions, which are related to the 
unpolarized parton densities; the latter are given by the MRST01 set [Hj], and the fragmentation functions by the 
KKP set [36j] . The valence u and d quark Sivers functions are the same as extracted from E704 data in Ref. [l[. 

Our results are shown in Figs. [JJ and [U and require some comments. 

• The thin, red, solid line in Fig. [TJ is the result of computing An using only the valence u and d Sivers functions, as 
extracted in Ref. [l|. 

• The cyan, dot-dashed curve in Fig. [TJ corresponds to the extreme case of the largest (in magnitude) gluon Sivers 
distribution, obtained by saturating the natural positivity bound, see Eq. ([J), 

A N f g/p ,(x,k ± ) = -2f g/p (x,k ± ). (8) 

The sea-quark Sivers functions are assumed to vanish. This leads to a SSA definitely in contradiction with the data. 
It also leads to a strong violation of the BSR. As expected, from PHENIX data a much more stringent constraint than 
the simple positivity bound for the (absolute value of the) gluon Sivers function can be obtained. The maximized (in 
magnitude) GSF in Eq. ([8]) has been chosen negative as the data hint at a possible small negative value of An- The 
opposite choice would lead to very similar conclusions, actually with an even stronger disagreement with data. 

• The thick, red, solid curve in Fig. [JJ has been obtained still assuming that there is no sea-quark Sivers contribution, 
and looking for a parameterization of A N f g / p i yielding values of An falling, approximately, within one-sigma deviation 
below the lowest pt data. This somehow corresponds to the largest (in magnitude) acceptable gluon Sivers function, 
in the sense that any larger GSF would cause the SSA to lie outside the error bars. 

The corresponding x-dependent part of the GSF, normalized to its positivity bound, \A N f g / p i (x)\/2f g / p (x), is 
shown as the red, solid curve in Fig. H as a function of x. PHENIX data on An clearly impose a stringent upper 
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bound on the magnitude of the gluon Sivers function in the region of small x, where gluons play a crucial role. The 
constraint is much less significant at larger x values, as it is natural: there cannot be any strong correlation between 
gluon distributions at large x (where they are anyway negligible), and a physical observable, like An at mid-rapidity, 
which is mainly sensitive to small x values. 

The GSF obtained here, the red, solid curve of Fig. [2 leads, within the x range covered by the data, to a strong 
violation of the BSR. This could be avoided by imposing an even smaller (in size) GSF, giving a strong role to the 
BSR. 

• As a last attempt, we have released the assumption of vanishing sea-quark Sivers distributions. With the aim of 
exploring how large a GSF can be, we have done that in an extreme scenario, i.e. assuming that all sea-quarks 
(u s , u, d s , d, s, s) have a non-vanishing positive Sivers function which saturates the positivity bound [that is, 
A N f q / p t (x, k±) = 2fq g / P (x, k±)}. Their contribution could then cancel the negative contribution to An of a possibly 
large GSF. 

Again, we look for the largest negative GSF which, together with a positive maximized sea-quark contribution, 
gives a SSA approximately lying within one-sigma standard deviation below the data points, represented again by 
the thick, red, solid line of Fig. [TJ This curve results now as the sum of the (maximized) sea and valence quark 
contribution (blue, dotted curve) and of the new GSF (green, dashed curve). 

The corresponding, normalized, new GSF is plotted as the green, dashed curve in Fig. O It is the largest gluon 
Sivers function compatible with PHENIX data: actually, it is strongly artificially enhanced by the extreme (or 
rather, unrealistic) assumption about the opposite contribution from saturated (and all summing up) sea-quark 
Sivers functions. Nevertheless, it still indicates a rather modest GSF in the important small x region. Notice that 
the opposite choice, a negative contribution from the sea-quark Sivers distributions balancing a positive one from the 
GSF, would lead to an even smaller GSF. 

Finally, we notice that, concerning the BSR, within the x range covered by the data, the (over)maximized sea-quark 
Sivers distributions give a positive contribution which strongly suppresses the negative contribution of the GSF, so 
that in this scenario the BSR is satisfied within a 10% level accuracy. 



III. COMMENTS AND CONCLUSIONS 



Unintegrated, TMD parton distribution and fragmentation functions may help in explaining several puzzling mea- 
surements for spin observables in inclusive particle production at high energy and moderately large px (for hadronic 
collisions) or Q 2 (for SIDIS). In particular, the quark Sivers distributions and the Collins fragmentation functions 
have recently raised a lot of interest. Azimuthal and single spin asymmetries measured by several experimental 
collaborations unambiguously indicate that these effects are sizeable, at least in some kinematical regions. 

A combined analysis of inclusive pion production in hadronic collisions and in SIDIS has given useful information 
on the quark Sivers distributions in the valence region. On the contrary, sea-quark and gluon Sivers functions are 
largely unknown. In this paper, we have performed a phenomenological analysis of the available data for the SSA 
An(p p ""A"), in particular regarding the recent PHENIX data obtained at mid-rapidity and moderate px values; 
this SSA is dominated by the Sivers effect and explores the low x region, resulting very sensitive to gluon contributions. 
To this end, we have adopted the theoretical approach of Refs. P,H, y|. 

A first result of our analysis is that all available data, including the PHENIX ones, are compatible with valence- 
like quark Sivers distributions and vanishing sea-quark and gluon contributions. The parameterizations required to 
reasonably reproduce the observed asymmetries also fulfill, within a 10% accuracy, the Burkardt sum rule. Since these 
parameterizations have a valence-like nature, the sum rule can be checked both over the fully integrated x range and 
in the limited x range effectively covered by data, with the same conclusions. 

The main issue of this paper was, however, that of discussing to what extent the available data can bind the 
magnitude of the gluon Sivers function, in particular in the small x region, where gluons play an important role. Our 
analysis shows that the PHENIX data on Am are presently the only ones that allow to reach quantitative conclusions 
on the magnitude of the GSF. The weakest upper bound (that is, the largest GSF), is obtained by balancing the 
gluon contribution to An with (over)maximized sea-quark contributions, opposite in sign. The resulting GSF is 
represented by the green, dashed line of Fig. O This scenario corresponds to a SSA lying approximately within one- 
sigma deviation below the data, as shown by the thick, red, solid line of Fig. [TJ Even in this extreme case, the bound 
obtained significantly reduces the magnitudes available from the simple positivity bound implicit in the definition of 
the Sivers function, at least in the x region where gluon contributions are relevant. 

Some comments are still in order: 
1) Following Refs. [3, Hi], 

we have assumed for the TMD functions a simple factorized expression, with a Gaussian 
k± dependence. The same (kj_) ~ 0.64 (GeV/c) 2 has been used for all partons (quarks and gluons). It has been 
suggested that gluons may have a larger (k±) than quarks [13] • We have found that, using smaller values of (kj_) q 
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and, e.g., (kj_) g ~ 2 (kj_) q would lead to a more stringent bound on the GSF. 

2) In our calculations we have used the MRST01-LO set for distribution functions and the KKP-LO set of fragmenta- 
tion functions. Since low- a; values are relevant in the kinematical configuration of RHIC experiments, one may wonder 
whether the choice of unpolarized PDF's and FF's affects the bounds. We have checked that there is indeed some 
residual dependence of this type. As an example, using the CTEQ6 PDF's [38|, which have a larger gluon component 
than the MRST01 set, would lead to a slightly more stringent bound. On the contrary, the use of the Kretzer FF set 
[3{|, which has a smaller gluon fragmentation component, would lead to a slightly less stringent bound. These small 
changes are well within the overall uncertainty of our results and do not alter our main conclusions. 

Future measurements covering the negative xp range, but at smaller cm. energies, e.g. at RHIC with yfs = 63 GeV 
or at the proposed PAX experiment at GSI, with y/s — 14 GeV, will certainly help in clarifying the size and relevance 
of the GSF. SSA measurements for inclusive photon production might also give useful information. A detailed analysis 
of the phenomenological interest of these processes is in progress [4(| . 
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Note added in proof 

Just after submission of our paper an interesting work appeared [4l| in which similar conclusions about the smallness 
of the gluon Sivers function are obtained from the smallness of the Sivers effect in SIDIS off a deuteron target. 
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PHENIX data on An- The red, solid curve corresponds to the case of vanishing sea-quark contribution; the green, dashed 
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